Exact form of the exponential correlation function in the glassy super-rough phase 
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We consider the random-phase sine-Gordon model in two dimensions. It describes two-dimensional 
elastic systems with random periodic disorder, such as pinned flux-line arrays, random field XY 
models, and surfaces of disordered crystals. The model exhibits a super-rough glass phase at l ow 
temperature T < T c with relative displacements growing with distance r as ([9(r) — (9(0)] 2 ) ~ 
A(t) ln 2 (r/a), where A(t) = 2r 2 - 2r 3 + C(r 4 ) near the transition and r = 1 - T/T c . We calculate 
all higher cumulants and show that they grow as {[0(r) - 0(0)] 2n ) c ~ [2(-l) n+1 (2n)!<(2n - l)r 2 + 
C(r 3 )] ln(r/a), n > 2, where C, is the Riemann zeta function. By summation we obtain the decay 
of the exponential correlation function as (e^M -9 ' )'} ~ (a/r)''' 9 ' exp ( — \A(q) ln 2 (r/a)) where 
r](q) and A(q) are obtained for arbitrary q < 1 to leading order in r. The anomalous exponent 
is r/(q) — cq 2 — r 2 q 2 [2'jE + ip{q) + ip(—q)\ in terms of the digamma function ip, where c is non- 
universal and 7_b is the Euler constant. The correlation function shows a resonant deep at q — 1, 
corresponding to fermion operators in the dual picture, which should be visible in Bragg scattering 
experiments. 

PACS numbers: 64.70.Q-,64.60.ae 



I. INTRODUCTION AND MAIN RESULTS 

The random-phase sine-Gordon model is the simplest 
model to describe the effect of quenched disorder on a 
periodic elastic system, the so-called random periodic 
class P"*' It models a host of experimental systems in 
presence of substrate impurities, such as charge-density 
waves, vortex latticesp^ random-field XY modelspl^HlI] 
smectic liquid crystals^ in situations where topological 
defects are absent, or can be ignored, allowing for an 
elastic description. In the simplest situation, the dis- 
placement fields from the perfect periodic position are 
encoded by a 2tt periodic one-component phase field 9(x). 
In a one-dimensional crystal of spacing a with elastic dis- 
placements u(x), the phase field is 9(x) = G u(x) where 
Go = 2ir/a is the smallest reciprocal lattice vector. In 
three dimensions, the random-phase sine- Gordon model 
was used to predict the existence of a Bragg-glass phase 
which is a glass pinned by (weak) quenched disorder, 
but which also retains topological or der (n o dislocations) 
and nearly perfect translational order^^ES] called quasi- 
order. Diffraction experiments^ measure the spatial cor- 
relations of the field e %q6 ^ x ' for q = 1. In that case the 
elastic description predicts a power-la w deca y of these 
correlations due to quenched disorder}- 8 1 13 1 17 1 character- 
istic of a quasi-ordered phase and leading to divergent 
Bragg peaks in the experiments. At stronger disorder, 
the presence of free topological defects is expected to 
lead to an exponential decay of these correlations with 
distance, as quasi-order is destroyed. 

Quasi-order usually arises when the phase field defor- 
mations 9(x) — 9(0) grow logarithmically with scale x, 
leading to power law scaling for the exponential field 
e iq6(x) _ rp p rODe deeper the properties of the phase field 
0(x) in a quasi-ordered phase it would be useful to pre- 
dict, and to measure, these correlations for q not neces- 



sarily an integer. One example with q — 1/2 would be a 
spin density wave, e.g., of XY symmetry S = Ae 1 ^ , sub- 
mitted to time-reversal invariant disorder, such as ran- 
dom anisotropy^, i.e., which couples to 2<j) = 6. Arbi- 
trary values of q would allow to probe the probability 
distribution of the phase deformations 8(x) — 9(0) and 
to characterize the multi-fractal properties of the expo- 
nential field. Here we restrict to the case of a d = 2 
dimensional periodic system, leaving the study of d > 2 
for an upcoming workPSl 

In two-dimensional periodic systems with quenched 
disorder thermal fluctuations play a more important role 
than in three-dimensional ones. As was discovered in 
the pioneering work of Cardy and Ostlund-^ they induce 
a phase transition at some critical temperature T c to a 
high temperature phase where disorder is irrelevant. For 
T < T c a glass phase exists in this m odel wh ich has been 
investigated in a number of works EH 4 I 19 I 20 I A very nice 
realization of this model in terms of crystal surfaces was 
described by Toner and DiVincenzo^L Since it allows in 
principle to measure the exponential correlation for any 
q, we now recall the basic phenomenology of surfaces. 
Note that another interesting realization of the Cardy- 
Ostlund model was obtained recen tly in the context of a 
smectic with surface disorder! 21 ! 22 ! 

Perfect crystals are characterized by an ideal lattice. 
At high temperatures, the thermal motion of the atoms 
overwhelms the lattice potential, and the crystal melts. 
In the present article we consider physical effects that oc- 
cur at significantly lower temperatures. Consider atoms 
at the surfaces of a crystal. They can more easily be dis- 
placed from their equilibrium positions by thermal fluc- 
tuations, since they reside at the boundary between the 
crystal and, usually, a fluid. They feel the periodic poten- 
tial created by the bulk of the crystal, as well as a more 
uniform potential from the fluid. At low temperatures, 



surface atoms are not displaced significantly from their 
equilibrium positions determined by the bulk potential, 
the surface is flat, and the atoms are ideally arranged. 
At higher temperatures, the thermal motion of surface 
atoms becomes significant, and they are randomly dis- 
placed, forming a rough fluctuating surface. The periodic 
potential is destroyed by thermal effects, a phenomenon 
known as the roughening transition! 23 ! I 

On the other hand, real crystals always experience 
some kind of disorder that tends to diminish the infinite 
correlation length of the translational order of a perfect 
crystal. In such situations, atoms at the surface do not 
longer experience a perfect periodic potential, but rather 
a disordered one created by the bulk. At high temper- 
atures the disordered potential is washed out by ther- 
mal fluctuations and thus unimportant for the shape of 
the surface: The surface is rough. On the contrary, at 
temperatures T below a critical temperature T c , surface 
atoms follow the disorder potential, thereby forming a 
surface that is even rougher. This phase transition is 
known as the super-roughening transition. 2 

Using surface-sensitive scattering experiments one can 
probe the crystal surface. One directly measures the dis- 
order and thermally averaged correlation function 



C(q,r) 



z iq{6(r)-e(0)]\ 



(1) 



where 9 denotes the two-dimensional height field of the 
surface, in units of Go = 2n/a. The quantity q denotes 
the component of the wave vector perpendicular to the 
surface. 

In the high-temperature rough phase, at T > T c , the 
correlation function C(q, r) decays at large distances as a 
power law C(q, r) ~ (a/r) r? ^\ where a is the lattice con- 
stant. The naive argument is that at high temperatures 
the fluctuations of the surface are effectively Gaussian, 
which in two dimensions always produce quasi-long-range 
order characterized by a power-law decay of C(q, r) and 
77(g) ~ q 2 . However, despite being an irrelevant operator 
at high temperatures, the lattice potential still plays an 
important role and produces in fact a nontrivial result 
for 77(g). In their seminal paper, Toner and DiVincenzo- 
calculated C{q, r) in the high-temperature rough phase, 
for T > T c . They obtained the decay as a superposi- 
tion of power laws C(q,r) ~ Ylm=i ^2m( a / r )'' m with 
exponents 

l lm (q)^(c-b)q 2 + b[(q- m ) 2 +m] . (2) 

At sufficiently large r the smallest power dominates, 
hence they concluded that the asymptotic exponent 77(g), 
for q > 0, takes the form 



77(g) = nun r] m (q) 

meN 



eg 2 + 6 [g] 2 + [g](l-2g) 



(3) 



where [g] is the integer part of q. The coefficients of 
Eq. (|3) are 

»-?. (4) 
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FIG. 1. The function 77(g) given in Eq. Q for b = 3, c = 3.2 is 
represented by the solid line. The functions 77 m (q) of Eq. H 
for m = 0, 1, . . . , 4 are given by the dashed lines. 



They are given in terms of the critical temperature T c 
and a non-universal off-diagonal disorder a' > defined 
in Sec. HH While the coefficient b is a universal function 
of the temperature, the coefficient c depends on a' and 
hence is non-universal. (It depends on the details of the 
model at short scales, which can renormalize a'). The 
form of Eqs. pi, pi has a simple interpretation in the 
picture of the Coulomb gas put forward by Cardy and 
Osthmd'3 (see also below). In that picture, perturba- 
tion theory to order 1p in the disorder is equivalent to 
inserting p pairs of replica ±1 charges. The vertex oper- 
ator e tq [ e ( r '~ 6 (°'' can itself be seen as inserting a charge 
+q at r and a charge —q at in the same replica. For 
q > I, it is then energetically advantageous to screen 
the vertex operator e*^ M ^'J with 2m replica charges 
from the disorder leading to Eqs. pi, pi 25 . As a conse- 
quence, 77(g) always grows with q, is quadratic for non- 
integer g, and continuous but with a cusp at integer val- 
ues of g. An interesting possibility, suggested by the form 
of the perturbation theory of Toner and DiVincenzo'3 
(i.e., the A-x m appear to have alternating signs) is that 
C{q,r) is an oscillating function, changing sign at inte- 
ger q. A simple toy model with such an oscillating be- 
havior is given in Appendix [Gj Finally, it is important 
to note that for a perfect crystal (i.e., in the presence 
of commensuration only and with no disorder), one has 
77(g) = b(q— [1/2 + g]) 2 , i.e., a periodic function of g with 
cusps at half integers and minima at all integers. 2 

Let us now consider low temperatures, T < T c . In this 
super-rough phase, the surface becomes even rougher 2 
producing a faster decay for C(g, r) as a function of the 
distance r. Although the general form for the decay of 
C(q, r) was correctly anticipated in Ref. 26 in the context 
of an A-componcnt extension of the random-phase sine- 
Gordon model, its precise expression, including the value 
of the exponents, has not yet been obtained. The aim of 
the present article is to fill this gap and determine the 
exact form for C(q,r) in the super-rough phase, close to 
the super-roughening transition. We will show that it 
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FIG. 2. Plot of the universal part of the anomalous exponent, 
rj(q), for < q < 1. 



takes the form 



/a\ r i(i) ( 1 „ \ 

G(g,r)~(") exp(^--.4(< Z )m 2 (r/a)J 



(5) 



characterized by the anomalous exponent 77(g) and an 
amplitude A(q). We obtained the exact dependence on 
q of both quantities, to leading order in t, the distance 
from the transition, defined as 



r = 1 - T/T c 



(6) 



For < q < 1, we find, up to terms of order C(t 3 ), 

v(q) = c q 2 + v(q), ( 7 ) 

V(q) = -T 2 q 2 [2 lE + TP(q) + iP(-q)] 1 (8) 

A(q) = 2q 2 T 2 . (9) 

Here i/j(q) is the digamma function. By rj{q) we denote 
the universal part of the exponent plotted in Fig. [2j it 
starts at 0(q 4 ) at small q; together with A(q) they are 
universal. The coefficient c in rj(q) is non- universal, and 
within our cutoff scheme given by 



= 2(1 -t)[1 + 2(1 -t>'] 



G(r 3 ), (10) 



where a' denotes a non-universal off-diagonal disorder 
defined below. 

Both the amplitude and the exponent are functions of 
the temperature T, with A = for T > T c ; r/(q) matches 
smoothly at T c to the expression ffib, i.e., rj(q) vanishes 
for T > T c . An initial guess for the amplitude may be 
obtained from the variance of the height fluctuations at 
two points, which we have recently calculated in Rcf. 19 
to two loop accuracy, 

{[9(r) - 6>(0)] 2 ) ~ A(t) \n 2 (r/a) + 0( ln(r/o)), (11) 

A(t) = 2t 2 -2t 3 + C(t 4 ). (12) 

This result was tested in a numerical simulation in 
Ref. 1201 using a dimer-model representation of the sine- 
Gordon Hamiltonian. If the displacement fluctuations 



A(q) 



Gaussian 



A{r)q 2 



(13) 



for all q. Interestingly, the present more detailed cal- 
culation confirms that property for q < 1. For q = 1, 
the amplitude A(q) of the correlation function C(q, r) 
changes abruptly to 



A(q = 1) = 6r 2 



(14) 



instead of A(q = 1~) = 2r 2 . Hence our results both for 
rj(q), and for A{q — 1), show deviations of the probabil- 
ity distribution of 9 from a Gaussian. Since the functions 
77(g) and A(q) are both increasing functions of q, the cor- 
relation function C(q, r) is a decreasing function of q for 
< q < 1. The amplitude A(q) jumps from 2t 2 at 
q = 1~ to 6t 2 at q = 1 resulting in a much faster decay. 
As a precursor of this effect, rj(q) diverges as q approaches 
unity. Such a resonance should be visible in Bragg scat- 
tering experiments, once the scattering wavevector or- 
thogonal to the surface matches Go = 2ir/a. 

For q > 1 and T < T c , one naively expects, both for 
rj(q) and A(q), additional resonances at wave vectors that 
are integer multiples of Go, and that the screening mech- 
anism which operates for T > T c in Eq. pi is also im- 
portant there. While a preliminary study indicates that 
this is the case, the detailed study of q > 1 is more com- 
plicated and deferred to a future publicationP3 

It is interesting to note that the case of integer q is 
relevant for the system of tw o-dimensional free fermions 
in a disordered potential ! 28 * 29 ! This model and the present 
random-phase sine-Gordon model are in correspondence 
via bosonization. More precisely, C(q — l,r) can be 
obtained as a four-fermion correlation function. We have 
performed that calculation, and found the result to agree 
with our expression ( 14) for q = 1. A more general study 
of Aq fermion correlator, that would allow us to study 
higher integer values of q in the correlation function, is 
in progress. The calculation and results will be presented 
elsewhere P3 

This article is organized as follows. In Sec. [n] we de- 
fine the model. In Sec. Mil wc introduce a formalism that 
enables us to produce a controlled expansion for the cor- 
relation function of interest. In Sec. II VI we evaluate this 
correlation function. This is followed, in Sec. |VJ by a 
discussion of the consequences, and conclusions. The de- 
tailed calculation of several involved integrals and some 
other technical details are presented in Appendices. 



II. MODEL AND THE PHASE DIAGRAM 

We consider the two-dimensional XY model for a real 
displacement field 9{x) € (—00, 00), without vortices, and 
in the presence of a random symmetry-breaking field. In 
the realization of the model as a fluctuating surface of 
the crystal, 9 describes the two-dimensional height of the 
surface. We have already considered the same model in 



an earlier publication^ here we repeat the necessary def- 
initions for the present study. The model is defined via 
its Hamiltonian 



H = d 2 x 



(V a 



h- V x t 



1 



(t* 



h.c. 



(15) 



where k is the elastic constant, a the lattice constant 
which provides a short-length-scale cutoff, and h(x) and 
£ (x) are quenched Gaussian random fields, the first one 
real and the other complex. Their nonzero correlations 
are given by 



h i (x)hJ{y)=T*—S ,3 6(x-y), (16) 

Z7T 

(17) 



Z7T 



where i, j € {1, 2} denote the components of h and T is 
the temperaturePB Note that the disorder h(x) must be 
introduced as it is generated by the symmetry-breaking 
field under coarse graining. We denote disorder averages 
by an overline. Depending on the context, x and y will be 
used cither to denote two-dimensional coordinates (as in 
the previous equations) or as their norms, i.e., x stands 
either for x or |x|. We emphasize here that the argument 
of the exponent in ( 15 ) should in principle be e lG ° e with 



Go being the smallest nonzero reciprocal lattice vector 
of the crystal normal to its surface under consideration. 
For simplicity we have set it to unity, thus measuring the 
displacement field in units of 1/G . 

We use the replica method to treat the disorderPH In- 
troducing the replicated fields 6 a , a — 1 . . . n, where by 
Greek indices a, j3, ... we denote replica indices, the repli- 
cated Hamiltonian reads 



"o +J1 1 ! 



H rcp - 
with the harmonic part 

■■J2[d 2 x{^5 a p[(V x 9 a ) 

-7~ Vcc# a • V x 9j3 >. 



(18) 



■"0 

T 



m 2 (6 a ) 



a/3 



(19) 



The mass m is introduced as an infrared cutoff. Wc 
perform calculations with finite m and study the limit 
m — > at the end. The system size is infinite throughout 
the paper. The anharmonic part reads 



H mp 



<] 



2ira 2 



y^ / d 2 xcos(9 a —6p). 

a{3 



(20) 



We start by c om puting the correlation function for the 
harmonic part (19), i.e., for g = 0. One obtains^ 



G a p(x) = (9 a (x)0ft(0)) = 8 afs G{x) + G (x), (21) 



where (. . .) denotes an average over thermal fluctuations, 
while at small distances \x\ <C mT 1 we obtain: 

G(x) = -(l-T)\n[c 2 m 2 (x 2 + a 2 )] , (22) 

G (x) = -2cr(l - r) 2 In [ec 2 m 2 (x 2 + a 2 )] . (23) 

Here we have introduced the ultraviolet regularization by 
the parameter a and c = e 7B /2 with 7^ being the Eulcr 
constant. The dimensionless parameter 



r = 1 - T/T c 



(24) 



measures the distance from the critical super-roughening 
temperature 



T r . = Ann. 



(25) 



The model studied here possesses an important sym- 
metry, the statistical tilt symmetry (STS), i.e., the non- 
linear part H[ cp is invariant under the change 9 a (x) — > 
9 a {x) + <f>(x) fo r an arbi trary function 4>{x). As discussed 
in many works B B I 32 I 33 I this implies that Gq(x) does not 
appear to any order in perturbation theory in g in the 
calculation of, e.g., the effective action.^ 

Let us summarize the one-loop re normalization group 
equations for the model ( 15 (j 1 ! 3 ! 4 ! 19 ! In terms of the scale 
£ := — In m they read 



dr 



= 0, 



dgii 
d£ 



2rg R - 2g R , 



d(TR 
dt 



1 



9%, 



(26) 

(27) 
(28) 



where the subscript r denotes the renormalized parame- 
ters that flow with the scale. Equation (26) is an exact 
result at all orders due to the above mentioned STS. We 
also see that the parameter a R does not enter any equa- 
tion (again due to STS), apart from being created by g R 
as one goes to larger scales. From Eq. (27) one reads off 
that the critical temperature is at r = 0. For T > T c we 
have scaling of gu to zero at large scales which denotes 
the rough phase. At low temperatures T < T c one finds 
a line of nonzero fixed points for g R that determines the 
super-rough phase. We see from the above scaling equa- 
tions that an grows unboundedly at all scales due to the 
non-zero fixed-point value for g R ; hence a R increases the 
two-point correlation function from a rough logarithmic 
behavior in the high-temperature phase to a super-rough 
square-logarithmic form for low temperatures. 2 In Fig. [3J 
we show the phase diagram for the model ( p~5J ) . The 
equations at two-loop order have been derived and stud- 
ied in Ref. [19j They do not quantitatively change the 
conclusions from the one-loop analysis, but give further 
insight into the two-point correlation function in the low- 
temperature phase. 
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FIG. 3. Phase diagram for the model (151. The high 



temperature rough phase is characterized by zero fixed point 
values of the anharmonic ga term, while the low-temperature 
super-rough phase has a line of fixed points for gn having fi- 
nite values. We are primarily interested to evaluate C(q, r) in 
the low temperature phase, for T < T c . 



III. GENERATOR OF CONNECTED 
CORRELATIONS 

The main goal of the present study is to compute the 
exponential correlation function 



C(q,r) 



s iq[6(r)-8(0)]} h 
-/g^P-yM-MO)] 



H"=p 



for the model (151 in the low-temperature phase. For 



actual calculations we use the second line of Eq. ( 29 1 that 



is obtained using replicas. By 7 we denote a particular 
replica index. 



In order to calculate the correlation function ( 29 ) it is 
useful to calculate the generator of connected correlations 
W(J). It is defined asP 



,W(J) 



V\ 



,-H( v )/T+J<p 



(30) 



where H(<p)/T = H ((p)/T+gV(ip) is the reduced Hamil- 
tonian. Here Hq is the quadratic part and V some pertur- 
bation. After performing the shift of the field as ip = x + 
GJ, one obtains H n (<p) — (pJ — H n (x)~H (GJ), where G 
is the propagator defined by H /T — ~((p\G \cp) where 
we use the compact notation. Equation ( 30 ) then be- 



,W(J) _ 



X>xexp 



flo(x) , H (GJ) 



T 



T 



-gV(GJ + X ) 



-Zq exp 



H (GJ) 

T 



(exp[-gV(GJ + X )}) x - (31) 



Using the cumulant expansion we obtain a perturbative 



expansion for W(J) that reads 



W(J)=\nZ 



<J 



H {GJ) 
T 



<(V{GJ + X )Y 



+ • > (v 2 (GJ + x)) X c +0(g 3 ). 



(32) 



In the following we calculate a perturbative expansion 
for the functional W(J) using the previous expression 
(32 1. Having evaluated W(J), one can immediately ob- 



tain all connected correlation functions by differentiat- 
ing the left- and right-hand side of its defining equation 
[equivalent to Eq. poll 



(33) 



W(J) =\nZ + In (cxp{Jip)) H 



an arbitrary number of times with respect to the source 
field J and setting J to zero at the end. For potentials 
that are even, i.e., when V(ip) = V(—tp), only terms with 
an even number of if fields have non-zero correlations. 



A. Expressions for W(J) 



Using the derived formula ( 32 1 and applying it to our 



(29) replicated Hamiltonian ( 18 ) we obtain 



W(J) = In Z Q + 



T 



W 1 + W 2 +0(g 3 ), (34) 



where the term Wj comes from the corresponding term 
of Eq. ( 32 ) proportional to gi . Either directly calculating 
or using the obtained results for the effective action from 
Ref. 1191 and the correspondence from Appendix [5] we ob- 
tain the final results for W\ and W%. We note that in 
models that have STS it was shown that two-replica part 
of the functional W and of the effecti ve act ion T are iden- 
tical up to the replacement 6 — > C/j 33 * 35 * see Appendix 

E 

The lowest order term reads 



W x 



5_ e -G(0) 



27TO 2 
X 



>J / cPxcos < 

a/3 



d 2 yG(x- V )[J a {y)-Jp(y)_ 



=0. 



(35) 



The last line is obtained in the limit we need, that is 



0, when the power law e 



-G(0) 



(c 2 m 2 a 2 y 



pensates the logarithmic divergence of G(x) under the 
cosine. The second order term W2 contains several dif- 
ferent contributions. In Appendix |A| we give a complete 
expression. Here we use only the final result taken in the 
limit m — > 0. It reads 



w 2 



8, 2 a^J 



d xd y 



(x — y) 2 + a 2 



2(l-r) 



d z z 



(l-r)ln 



{y-z) 2 + a 2 
(x — z) 2 + a 2 



[J a {z) - Mz)} 



(36) 



W = ln 



-D- 



W 2 = -g 




-m—o- 



FIG. 4. A graphical representation of the non-vanishing 
contributions to the perturbative expansion. The white 
and black squares stand for the observable e lq6l ^ r ' (with 
charge +q) and e~ lq ~* °' (with charge —q) respectively. The 
replica indices are indicated by a solid line. The interaction 
Enifl 6 '^"' 1 ' - ^' 1 " ls coded with a white and black circle, 
joined by a dashed line. There has to be as many black as 
white points and circles in each given replica (charge neu- 
trality). Charges qi at x\ and qi at Xi interact with a mul- 
tiplicative factor (Boltzmann weight) of \x\ — X2\ 29192T ; 
they "attract" for opposite signs, leading to a divergent con- 
tribution as their distance goes to zero. This allows for an 
intuitive visualization of the corresponding divergences in the 
perturbation expansion. 



B. Special source field 

Using a source field of the form 

Ja{%) — iq [8(x — r) — 8{x)] S ai , 



(37) 



in Eq. (33), we obtain a relation between the correlation 
function (29) and the functional W(J a ). It reads 



M0-,{r)-0-,m 



1 



,W(J a 



H vap 



(38) 



The left-hand side of Eq. ( 38 ) is the correlation function 



C(q, r) , while the right-hand side of it should be evalu- 
ated. Here we perform a perturbative calculation in the 



anharmonic coupling [defined in Eq. (20)] of W(J a ), us- 
ing the expansion p4| 



For the source field (37), the lowest order term in the 



expansion of W{ J a ) becomes 

K P (GJ) _. 

T 

a/3 

= -q 2 {\ - t) [1 + 2(1 -T)a] In 



W 



£ / d?xd?yJ a {x) Ga ^l y) J fi Un 



r 2 + a 2 



(39) 



A graphical interpretation is given in Fig. [4] 

The second nonzero term of W{J a ) is obtained by em- 
ploying the source field (37) in Eq. p6|, and by mak- 



ing use of the summation formula over replica indices 



Z)™ «=! COs[A(8 ai - 8 fa)] =n 2 +2(1 -n)(l -cos A). In 
the limit n — > 0, we obtain 

W 2 = -g 2 {^y T I(r,a,q,r), (40) 

where we have defined the following integral: 



I(r,a,q,r) 



4tt 2 



d xd y- 



1 

[a 2 + (x - y) 2 ] 



2(1- 



[(x- 


r) 2 + a 2 y 2 + a 2 ' 


L(y- 


r) 2 + a 2 x 2 + a 2 



q(l-T) 



(41) 

It is an even function of q, and as a function of r and a 
depends only on r/a, provided it is convergent. 

Finally, using Eqs. @, @, @, and @, as well 
as the result Zq/Z = 1 valid in the replica limit n — > 
0, we obtain the final (yet unevaluated) result for the 
correlation function ( 29 ) , 



C{q,r) =exp -q 2 {\ - r)[l + 2(1 - r)a] In 



x exp 



-g 2 (^) 4T I(r,a,q,r)\ exp (G( 5 3 )) 



(42) 



Equations (41 )-(42 1 are the starting point of the eval- 



uation of the correlation function ( 29 1 . One should read 



them as the result of perturbation theory to leading or- 
der in the bare coupling g. However, if we re-express g in 
terms of the renormalized one gn 1 we obtain the result to 
leading order in r, the distance from the transition. They 
are connected by the relatiorPS g = g R [l + O(r)]. In the 
massless limit, equivalently, at large distances, r ^$> a, 
from the renormalization group equation ( |27[ ) the renor- 
malized coupling reaches its fixed point value t + 0(t 2 ), 
hence we can replace 



0{t 2 ) 



(43) 



in the evaluation of the large distance behavior. The be- 
havior of gn at intermediate scales before reaching the 
fixed point causes some non-universal behavior. It is 
easy to see, however, that it is confined to subdominant 
terms, such as the part of ln(r/a) that is proportional to 
q 2 . As a consequence of flowing of gn toward the fixed 
point, the off-diagonal disorder becomes changed at large 
scales. Therefore by a' we denote the non-universal pa- 
rameter that characterizes the off-diagonal disorder. We 
note that, to lowest order in g, we have a' = a. In order 



to emphasize the difference between the bare parameter 



a, defined in Eq. (16 1, and the effective one at large scales 
a', we keep a' in final formulas, e.g., in Eq. (10 1. 

Thus, to leading order, we only need the integral 
I(r, a, q, t = 0), that is a function of the ratio r/a and q 
only. It is clear from the definition (41 1 that for r = 



this integral is convergent for any q: It is ultraviolet con- 
vergent thanks to the cutoff a and infrared convergent 
thanks to the substraction of 1. There are however ul- 
traviolet divergences when a/r —> 0, which lead to a loga- 
rithmic behavior, and which will be analyzed in the next 
section. Here we want to point out that while for \q\ < 1 
they come from the first power law factor, there is an 
additional ultraviolet divergence arising from the second 
factor (containing r) when q = 1. Hence we will mainly 
restrict to q < 1, and discuss separately the cases |g| < 1 
and q = 1. It makes physical sense that at q = 1 (and 
in general any integer q) the correlation function changes 
non-smoothly, as was already the case in high tempera- 
ture phase [see the results (J3J) of Toner and DiVincenzo'21] . 
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FIG. 5. Numerical evaluation of Eq. (451 for j — 2, 3. At 



large r, we observe that the data belong to a straight line on 
a log-linear plot. The numerical coefficients [2£(3) and 2£(5) 
for j = 2 and j = 3, respectively] in front of the logarithms 
in the fitted straight lines are analytically obtained by the 
dimensional method, see the main text. 



IV. EVALUATION OF THE CORRELATION 
FUNCTION 

In this section we evaluate the correlation function 



where we have defined the family of integrals 



C(q, r) defined by Eq. ( 29 ) , using its form ( 42 ) , for gen- 
eral values of the parameter q in the region < q < 1. 
To achieve this, our task is to calculate the integral of 



Eq. (41 1. Exact evaluation of that integral is rather diffi- 



cult. However, we only need the most dominant term for 



large r/a, when the expression of Eq. (41 ) is expected to 
take a universal form. 

We calculate the integral using two different methods. 
The first one, which we term the finite-a method, works 
directly at r = 0, but with a non-zero ultraviolet cut- 
off a > 0. However we are not able to obtain all the 
results using this finite-a method. Hence we also use di- 
mensional regularization, which we term the dimensional 
method; it is quite powerful, although a little delicate in 
its interpretation. In that method one works directly at 
a = 0, but keeps r > 0, which renders the integral con- 
vergent. From the poles in 1/r one extracts the desired 
information. A careful comparison between the obtained 
results is performed at the end. 

A. Finite-a method 

We start the evaluation by setting r = and keep- 
ing a finite in Eq. (41). This is justified as we have 
I(r, a, q, t) — I(r, a, qjT) + 0(t), so to lowest order in the 
distance from the transition we only need I(r, a, g, 0), as 
discussed above. 

For small g, we expand it into a Taylor series as 



I(r, a, g ,0)=Vft ir/a) 



(44) 



Ij(r) 



1 



4tt2(2j)! 



d 2 xd 2 y 



1 



[i + (x-y) 2 r 



x In 



2j 



(x-r) 2 + l y 2 



(y-r) 2 + l x 2 + l 



(45) 



The expansion ( 44 ) contains only even powers of g, a con- 
sequence of the parity I(r, a, g, 0) = I(r, a, — g, 0). Rather 
involved evaluation of the lowest order term I\(r) of 
Eq. (44) is presented in Appendix [d] At small g, the 
result reads 

I(r, a, q, 0) - g 2 In 2 (£j + g 2 In (~) + r(g), (46) 

where the even function r(q) at small g starts with a term 
proportional to g 4 and contains all the desired higher- 
order terms proportional to g 2 - 7 , j > 2. 



To determinate r(q) defined in Eq. (46) we begin by 
a numerical evaluation. In Fig. [5] we show the results 
for Ij(r) for j = 2, 3 that determine the two lowest-order 
terms of r(q) at small q, see Eqs. Ell) and (46). We 
notice that the data points for 12(f) and ^(r) appear 
to be on a straight line on a log-linear plot, which shows 
that these two functions are well described by a logarithm 
ln(r) with j-dependent prefactors [2£(3) and 2£(5) in the 
present case]. This leaves us with a hint that all Ij(r) for 
j > 1 might asymptotically have a logarithmic behavior 
in the limit r>l. That would determine the anomalous 
exponent for the C(q, r) correlation function. 

We gain further knowledge about I(r, a, g, 0) by con- 
sidering the special case q = 1. The calculation presented 
in Appendix [D] reveals the result 

I(r,a, 1,0) = 3 In 2 (£) - 2 In (£) +O (j^j , (47) 



showing that the coefficient q 2 of the leading squared 



logarithmic term from the small-g expansion of Eq. ( 46 ) 
changes at q = 1. As discussed above this is expected 
from physical considerations. To confirm this result for 
q = 1 we have performed a calculation in the fermionic 
version of the model that allows us to treat integer q 
only. It confirms our result for q = 1. The calculation 
and results will be presented elsewhere.^ Note that at 
q = 1 the prefactor in front of ln(r/a) term in Eq. (47) 



takes this value for the particular cutoff procedure we 
use, but is not expected to be universal in general. 

Motivated by the above analysis, we introduce the 
anomalous exponent 77(g) and the amplitude A(q) and as- 
sume a general form for the correlation function at r >a 
given by Eq. |5|, where 



77(9) = 2q 2 (l - r)[l + 2(1 - r)a] + g%(q) + 0(g 3 ), 



A(q) = a(q)g 2 + 0(g 3 ). 



(48) 
(49) 



For the unknown coefficients in the previous two expres- 
sions, so far we have established the following results: 



a{q) = 




(50) 
(51) 



B. Dimensional method 



To evaluate Eq. (41 ) using the dimensional method, we 



take the limit a — > of the integrand, i.e., we consider 

I(r,a,q,T)=I(r,0,q,T)+O(a 2 ) 

=7(1,0, <z,T) + 0(a 2 ), (52) 

where the second line is a trivial consequence of rescaling 
the coordinates by r, provided the integral is convergent. 
For t > the integral l(l,0,q,r) is ultraviolet conver- 
gent as long as \q\(l — r) < 1; it remains infrared conver- 
gent if t is not too large (r < 1/2). We thus first discuss 
I q I < 1. From the poles of the evaluated expression in 
the limit r — > we will infer the behavior of the integral 
at nonzero a, as is shown below. We rewrite 



1(1,0,8,7-) 



1 

47r2 



J(q(l 



where we defined 



J( P ,t) :=FP 



d x d 



y|x- 



|4(r-l) 



>r), 



|*-e||y| 

|y-e||x| 



(53) 



■2l> 



(54) 

Here e is an arbitrary unit-vector, and FP means "finite 
part" in the sense of dimensional regularization (in some 
domain of q this finite part is achieved by the subtrac- 
tion 1 in Eq. ( 41 ) , however it can be given a more gen- 



parameters q, t). The evaluation of the ensuing 2x2- 
dimensional integral is complicated and the details are 
presented in Appendix [C] Let us recall here the main 
idea, which goes back to Dotsenko and Fateev^E^l The 
integral (54) can be thought of as an integral over the 
complex plane, both for x and y. The integral over, say 
x over the complex plane can be decomposed into two 
real contour-integrals, over x and its complex-conjugate 
x. Noting that (54) can be written in the form 



d 2 xd 2 y(l - x) p y p x 2r - l {l - j/) 2(t_1) (1 - xy)~ p 

x (1 - i)Yi 2T_1 (l - 2/) 2(T ~ 1} (i - xy)~ 



(55) 



it is suggestive that deforming the contour integrals over 
x and x to lie on the real axis, the resulting integral will 
take the form 



J(p,T) = J2 B i J i J i- 



The BiS are phase-factors - typically deforming a con- 
tour around a branch cut gives a sine-function of 27r 
times the power at the branch cut. The Ji are second- 
generation hypcrgcomctric functions 3-F2, e.g., the first 
line of Eq. ( 55 ) will naturally lead to 



J t = [ dx [ dy(l-x)Py p x 2T - 1 (l-y) 2 ^- 1 \l-xy)-P. 
Jo Jo 

Note that the integral above is restricted, both for x 
and y, to the interval [0,1]; but the domains [— 00, 0] 
and [1, 00] also contribute; the latter can then be trans- 
formed back to the interval [0, 1], leading to integrals of 
the same form, but with different coefficients. This form 
can be recognized in the exact result given in Appendix 
[C) formula (C3|. 

The final result for < q < 1 and small r reads 



/(l,0, g ,T) 



g 

~8t 2 



7 2 [1 



27£ " 1>(q) ~ 1>(-q)] 



At 



O(r ) 



(56) 



where ip(x) = V (x)/Y(x) is the digamma function. The 
coefficients of the poles in Eq. ( 56 1 contain all informa- 



tion needed for the determination of the anomalous ex- 
ponent ( 48 ) and the amplitude ( 49 ) . The precise method 



eral meaning in terms of analytical continuation in the 



to extract the information from these poles is given in the 
Appendix |E 1| Here we give a more intuitive derivation. 
Starting from Eq. (40 ), which is contained in the cor- 
relation function ( |42[), w e expand (r/a) 4T I(r,a,q,T) at 
small t. Using Eqs. (|52[) and (56), for < q < 1 we find 
the result 

^ I(r,a,q,r) = q 2 \n 2 Q + Vg (q)ln(^j + 0(r), 

(57) 
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where the prefactor of the logarithmic term, i.e., the non- 
trivial part of the anomalous exponent, reads 

r) g (q)=q 2 [l-2 lE -^(q)-i>(-q)}. (58) 



Equation (57) could also be understood as the final ex- 



pression for I(r, a,q,T = 0) obtained via the dimensional 
method. At small q, one can notice the agreement be- 
tween the two results p6| and (57) that are obtained 
using quite different methods. 



The expression ( 57 1 should be taken with special atten- 
tion. The translation of the dimensional-method result 



( 56 1 by performing a naive expansion at small r of the 



left hand side of Eq. (57) would contain extra terms 



q 

8r 2 



1 

4t 



2q 2 In - 
a 



■Vg(q) 



(59) 



on the right hand side. These two terms are divergent 
for r — > 0. The terms (59) formally appear only because 



we performed a small a expansion in Eq. (52), i.e., we set 



a = in the integrand. The left-hand side of Eq. (57) 
should be finite at r = 0, and thus the right-hand one as 
well. Therefore, the extra terms (59) should be omitted 



when one relates the dimensional result ( 56 1 to the ex- 
pression for I(r,a,q,r = 0). We again emphasize that a 
more rigorous derivation is given in Appendix |E 1[ while 
Appendix |E 2| we give an illustrative example. 



V. FINAL RESULTS AND CONCLUSION 



the correlation function 



Using Eqs. (|42|) and (|57|) we obtain the final result for 

(60) 



C{q,r)=(e^ e ( r )~m]) H , 



at q < 1, for the model (15). It takes the form 
'a\ r i(q) 



C(q,r) 



cxp 



- l -A{ q )\n\r/a) 



where the amplitude is given by 



J2g 2 T 2 + 0(r 3 ), q<\ 
A[q) ~\6r 2 + 0(T 3 ) : q = l 



(61) 



(62) 



Here and below we use the replacement g — > r as dis- 
cussed in Eq. (43), i.e., the fixed point of the renormal- 
ization group. 



The anomalous exponent of Eq. (61 1 reads 



Here 7^; is the Euler constant. This produces the result 
([7]) displayed in the Introduction where we have sepa- 
rated the universal part, 77(g), of 77(g) that starts at order 
q , from the part proportional to q 2 that is non- universal. 
The anomalous exponent 77(g) rapidly increases as q 
increases from zero to one. Thus, C(q\,r) > C(q2,r) 
for < qi < (72 < lj meaning the correlation func- 
tion decreases as q increases. When q approaches unity, 
77(g) ~ 1/(1 — g) becomes very large and makes the cor- 
relation function C(q, r) decay much faster, see Eq. ( |6"T| ). 
This is a precursor of the more drastic effect which hap- 
pens at q = 1 where the amplitude A{q) jumps (from 
2g 2 r 2 at g — >• 1~ to 6t 2 at g = 1). These dips in the 
correlation function C(g, r) are a remarkable and unique 
feature of the super-rough phase. Note that near q = 1 
there is a growing length scale 



u 



aexp 



1 



i-g 



(65) 



belo w which the 77(g) ln(r/a) term in the exponential in 
(61) is larger than the asymptotic In (r/a) term. 
Our result (J64J) can also be expanded at small q as 



r, g (q) = q 2 + 2C(3)g 4 + 2C(5)g 6 + 0(q s ). (66) 

This expansion precisely determines the prefactors for 
the family of integrals of Eq. (45). For j > 2 we find 



Ij(r) = 2((2j - 1) ln(r-) +O(r ), which explains the data 
of Fig. [5] This provides a confirmation that we have 
correctly extracted the amplitudes from the dimensional 
method (at least for small q). 



Our result for the correlation function (60) enables 



us to calculate the leading large-distance behavior of 
all higher powers of the connected correlation func- 
tions in the super-rough phase, i.e., for T < T c 
(see Fig. [3]). Using the cumulant expansion formula 
ln(expL4]) = YjjL\{-A?)clfi> after expanding (the loga- 
rithm of) Eqs. (60) and (61) at small g, one obtains 



[0(r) - 9(0)} 



2j-l 



0, 



-l) j 



(2j)! 



?(r)-0(O)r 



-2t 2 ((2j - 1) In 



-Oir^r 



for j ' > 1. On the contrary, for j = 1 one finds 



([8(r)-6(p)Y 



2r 2 ln 2 



■O In 



(67) 



(68) 



(69) 



77(g) = 2g (1 - r)[l + 2(1 - r)a'] + r Vg (q) + 0(r ) Expression p| is the well-known result™ for the 



(63) 

where its nontrivial part determined by the anharmonic 
coupling reads 



model (15). This is yet another way of obtaining the 



Vg(q) 



q 2 [l~2 lE -iiq)-^(-q% q<\ 



(64) 



result that produced some controversies in the past, as 
discussed in Ref. 01 Confirming a recently obtained 
correction's to the prefactor of the squared logarithm 
in Eq. ( 69 ) using the present method woul d r equire ex- 
plicit evaluation of the 0(g 3 ) term in Eq. (42) that is a 



formidable task beyond the scope of the present study. 
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Our main results directly apply to some physical sys- 
tems, in particular to surfaces of crystals with quenched 
bulk disordei^ or to a vortex lattice confined to a 
plane.' 5 In particular, the structure factor S(q, k) = 
J d 2 r C(q, r) exp(ik r) in the super-rough phase at k = 
is analytic and finite. 2 However, we predict that it has 
sudden dips once the wave-vector q becomes an integer 
multiple of the reciprocal lattice vector normal to the sur- 
face Go (set to unity here) of the bulk crystal; we expect 
it to have characteristic jumps in the amplitude and the 
anomalous exponent not only at q = 1 but rather at any 
integer q. 

To summarize, the main results of the present study 
are given by Eqs. (60)-(64), or equivalently Q. They 



describe the behavior of the exponential correlation func- 
tion in the super-rough phase, defined by Eq. ( |60| . We 
found a characteristic universal jump of the amplitude 
of the squared logarithmic term, at q = 1, which occurs 
in conjunction with the drop of the anomalous exponent. 
The value of the amplitude at q = 1 is in agreement with 
the result in the equivalent fermionic disordered model 
and will be published elsewherepS It would be very in- 
teresting to perform a numerical simulation, along the 
lines of Ref. HUJ to test the results of the present work. 
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Appendix A: W2 term 

In this Appendix we provide the complete expression 
for the second-order term W% of the generator of con- 
nected correlations, see Eq. (34 1. It reads 



W 2 



\{^f*-* am Yj:j^y 



where 



s=± a/37 ' 

■ < -A 1 (x -y)S ai cos [(GJ a /3)(x)] 
+ [A s 2 {x-y)5 ai +A\{x-y)] 
x cos [{GJ a p){x) - s(GJ 70 )(y)] 



At(x) = 4(2 - e~ G{x) - e G ^), 

A\(x) = 1 + e 2sG ^ - 2e sG(x \ 

A\{x) = 2(e sG{x) - 1). 



(A2) 
(A3) 
(A4) 



Here we introduced the abbreviation 

(GJ a0 )(x) = Jd 2 yG(x - y) [J a {y) - J p (y)} . (A5) 

In the limit of zero mass, m — > 0, the only term that 
survives in Eq. ( Al I comes from the A 2 term and reads 



W 2 = -£-Y [ d 2 xd 2 ye 2G ^-^- 2G ^ 
8lra ^ J 



xcosi / d 2 z [G(x - z) - G(y - z)] \J a {z) - Jp(z)\ > . 

(A6) 

The difference of propagators in the previous equation is 
well behaved in the limit m —$■ 0. Using the propagator 
( 22 ) , we eventually obtain the expression ( 36 ) of the main 



text. 



Appendix B: Connection between the effective 
action T and W 



The final formula (32) resembles the expression (A17) 



of Ref. [19] for the effective action, provided in the latter 
one omits all the terms that explicitly contain the prop- 
agator G and one multiplies all the terms in Eq. (A17) 
by a factor —1 apart from the harmonic part Sq (corre- 
sponding to Ho/T in the present work). We note that 
in general, in models that have STS, two-replica part of 
the functional W and T are identical up to the replace- 
ment of the form 9 —¥ GJ and a global minus sign, see 
Eqs. (89) and (90) of Ref. M 

For the model studied in the present work, we calcu- 
lated the effective action T in Ref.UHl the term T 2 propor- 
tional to g 2 is given by Eq. (25) there. We notice that in 
Ref. EH we used restricted summation over replica indices 
(all replica indices different), while in the present work 
we do not impose such restriction. As a consequence, if 
one wants to use the theo rem of similarity between the 
functionals W and r j 33 l 35 l one should transform the ex- 
pression for r to contain summations over unrestricted 
replica indices only, before applying the theorem. On 
the contrary, if one directly uses the expression for F2, 
Eq. (25) of Ref. [19j its two rep lica part corresponds to 



the two replica part of Eq. ( Al ), provided one omits the 
linear term —pG(x) in F2 and performs the replacement 
9 — > GJ, or written more explicitly: 



J, (Al) 0*(x)->J2 JcPyGcyix-vWy), 



(Bl) 

9 a (x) - 9 (x) -> / d 2 yG(x - y) [J a (y) - My)} 

= (GJ af> )(x). (B2) 

In the last expression we explicitly used the fact that 
G a p{x) = S a pG(x) + Gq(x), so that only the diagonal 
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cr-independent part survives. This is one of the manifes- 



tations of STS where we explicitly see through ( 32 1 that 



all corrections to the correlation functions that come due 
to the anharmonic term (proportional to g) are a inde- 
pendent. 



Appendix C: Integral I(r,a = 0,q,r) 

In this Appendix we analyze and for r > calculate 
the integral 7(1,0, q,r) defined in Eqs. (41) and (52); it 
can be written as 



I(l,0,q,r) = ~J(q(l-T),T), (CI) 



where the integral J(p,r) was defined in Eq. (54) in the 
main text. 



The integral of the type ([54]) was studied in Refs. [38] 
and I3HI in the context of correlation functions in a two- 
dimensional Ising model. It can be expressed in terms of 
the hypergeometric functions, defined as 



pFq( a l 



,a p ; b\,..., b q ; z) 



(oi)fc . . . (Op)* z k 



^ (&i) fe • • • (Mfc k\ 



(C2) 



where (a)k = T(a + k)/T(a), while p and q are positive in- 
tegers. In the following we use results from Ref. 39 (which 
are equivalent to those of Ref. [35] once one employs some 



identities between the hypergeometric functions). It is 
important to note that the integral J(p, r) diverges for 
p = 1 (or more generally integer p) but that the formula 
below provides an analytic continuation in p away from 
the poles at integer p. In turns, it gives an analytical 
continuation for the original integral, which we will call 
7 CO nt(l, 0, 9, T )- We now need to discuss how this contin- 
uation is relevant for the physical problem studied here, 
and online the content of the next two subsections. 

Our main objective, as discussed in the main text, is 
to extract from the poles in 1/r 2 and 1/V of this inte- 
gral the logarithmic behaviors at large distance r. How 
to do that is explained in Appendix |E 1[ In addition we 
will restrict our study to \q\ < 1. As was discussed in 
the Introduction, the divergences that appears at q = 1 
have important physical consequences and are related to 
the screening mechanism unveiled by Toner and DiVin- 
cenzo^. To make progress for q > 1 one needs to reex- 
amine the entire perturbation expansion, which we defer 
to a future work. The considerations directly useful for 
the present work are presented in the first part of this 
Appendix. 

Still, it is not devoid of interest to also study (i) the 
analytical continuation I C ont(l, 0, g, r) for arbitrary q (ii) 
the finite parts of the integral, even if it will not contain 
the universal information sought here. These considera- 
tions are reported in the second part of this Appendix. 
1. Study useful for \q\ < 1 and the divergent parts 



Using the results from Ref. 39, for the integral (54) one 
obtains 



J(p,r) 



sin 2 (np) 



2 sin (np + 2-kt) 



{ [cos(2?rp) + cos(2?rp + Attt) - 2] B(l -p,l + p) 2 B{\ -p,l- 2r) 



x 3*2(1 - p, 1 - 2r, 2 - 2r; 2, 2 - p - 2r; l) 2 

+ [cos(4ttt) - 1] B{1 + P ,2t- l) 2 B(l + p, 2t) 2 3 F 2 ( P , 1 + p, 2r;p + 2r, 1 + p + 2r; l) 2 
- 4 sin(Trp) sin(27rr)B(l - p, l+p)B(l - p, 1 - 2r)B(l + p, 2r - 1)B(1 + p, 2r) 
x 3*2(1 - p, 1 - 2r, 2 - 2r; 2,2 -p - 2r; l) 3 F 2 (p, 1 + p, 2r;p + 2r, 1 + p + 2r; 1)}, 



(C3) 



where B(x,y) is the beta function. We need a simplified form of Eq. (C3) valid at small values of r. 
The excess of the hypergeometric function 3*2(0, b, c; d, e; z) is defined as d + e— a — b — c and for 



has to be positive in order for 3*2(0, &, c; d, e; z = 1) to be defined by its series (C2 ). A naive expansion 
3 F 2 (l-p,l-2r,2-2r;2,2-p-2r;l) = 3*2(l-p,l,2;2,2-p;l) + 0(r) 



1 the excess 



(C4) 



of one of the hypergeometric functions of Eq. ( C3 ) is not possible. Therefore, although the hypergeometric functions 



of Eq. (C3 1 are well defined, since they have positive excess 2r, we cannot perform a Taylor expansion around r = 0. 



In the following we use particular relations between the hypergeometric functions in order to transform Eq. ( C3 ) 



into a s uita ble form that can be Taylor expanded at small r. After transforming 3*2(1— p, 1 — 2r, 2 — 2r; 2, 2— p — 2r; 1) 
of Eq. (C3| using the one-term transformation rule^ valid for a, s > 



3*2(0, b,c;d,e;l) 



r(d)r(e)r(s) 

r(a)T(b + s)r{c + s) 



3F 2 (d — a,e — a,s;b + s, c + s; 1), 



= d + e — a — b — c, 



(C5) 
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one obtains the following result valid at p < 1, 
J (p. t) 7r 2 p 2 [cos(27rp + 4"7tt) + cos(27rp) — 2] 



4tt 2 



8sin 2 (27rr) sin 2 (7rp + 2ttt)) 
ttV (2t - l)(p + 2r)r(2r - l) 2 



+ 



2 sin 2 (7rp + 2irr))T{\ - p) 2 T(p + 2r + l) 2 
p 2 (l - 2t) 2 sin 2 (2^T)r(l +p) 2 r(2r - l) 4 r(l - p - 2t) 2 
47r 2 r(l-p) 2 r(l+p + 2r) 2 



3F 2 (l+p,l-2r,2r;l,2;l)^ 

3 F 2 (p, 1+p,2t;p + 2t,1+ P + 2r; 1) 3 F 2 (1 + p, 1 - 2r, 2r; 1, 2; 1) 



,F 2 0, 1 + P, 2t; P + 2r, 1 + p + 2r; if 



(C6) 



The expression (C6) has a positive excess in the hypergeometric functions under the two conditions r > and p < 1. 



The main advantage however of Eq. (C6| is that it is more convenient for Taylor expansion around r = 0. We need 



to expand the two hypergeometric functions of Eq. (C6). The first one, 3^2(1 +P, 1 — 2r, 2r; 1, 2; 1), can be expanded 

rp< 1. One obtai 

1 - ie - i>(—p) + 



safely, since it is well-defined in the expansion limit r — > for p < 1. One obtains 

1 



3 F 2 (l+p,l-2r,2r;l,2;l) = l + 2 



P 



t + 0(t 2 ). 



(C7) 



We employed the Gauss relation 2-^ 1(0, 6; c; 1) = r(c)r(c — a — 6)/[r(c— a)r(c — 6)] and used the fact that 3F2 re duce s 
to 2-F1 for r = 0, see the definition (C2). However, the expansion of the second hypergeometric function of Eq. (C6|, 
3F2O, 1 + p, 2r;p + 2t, 1 + p + 2r; 1), can not be done directly, since for t = it has zero excess. To overcome this 
problem, we use the following relation^] 



3 F 2 (p, 1+ Pi 2t; P + 2t,1+p + 2r; 1) = (1 ± P)(1 + 2t) 3 F 2 0, 1 + p, 2r;p + 2r, 2 + p + 2r; 1) 



1 + p + 2t 
P(l-2r) 

P + 2t 3 



-F2O. 1 + P, 2r; 1 + p + 2r, 1 + p + 2r; 1) 



(C8) 



which connects 3^0, 1 +p, 2r;p + It, I + p + 2r; 1) with two other hypergeometric functions that have an excess of 
1 + t and hence are well defined for r = 0. One then obtains for the expansion 



3 F 2 (p, 1 + p, 2r; p + 2t, 1 + p + 2r; 1) = 2 + 2 



2 7£ + 2V>0) + 



+ 0(r 2 ). 



Combining Eqs. (C7), (C9), and (IC6|, one obtains the expanded form of Eq. (C3|, 

,2 



^ =^ + £ P - 2 7£ - *(p) " #-*)] + ©(A 



(C9) 



(CIO) 



Finally, using Eqs. (CI) and (CIO), wc obtain the result reported in the text (56) 



It is interesting to notice that Eq. (C3| satisfies the relation J(p,t) — J(p, Tj2 — r) when r > and 1/2 — r > 0, 
i.e. for < t < 1/2. This significantly simplifies the Taylor expansion, since Eq. (C3) could easily and directly be 



expanded around r = 1/2, because the excess 2r of the hypergeometric functions is then always positive and hence 
they are well behaved. Performing the calculation in this way one obtains 



*? = WJ ^ + j^ P- - 2 ^ - *0» - *H»1 + ««»/» - '>">. 



(Cll) 



which is in agreement with the result (C10). This observation suggests that a more symmetric form of the expansion 
reads 



P 



P 



JQ.t) 

4tt 2 32r 2 (l/2-r) 2 t 8t(1/2-t) 



-2 7£ - </>0) - ^(-P)] + O(r (l/2 - r)°). 



(C12) 



2. Further study of the analytical continuation 



We should notice that although Eq. (CIO) is formally obtained from Eq. (C6) that is valid for p < 1, Eq. (CIO) is 



an even function of p [as well as the starting ones of Eqs. (54) and (C3|] and t 
the parity transformation. Hence we can define the integra 



rerefore Eq. (CIO) holds at any p by 



analytically continued to the whole complex plane: 



-front (1,0, q,r) 



q 

8r 2 



q 2 [l-2 lE -^(q)-jj(-q)} 
At 



+ 0(t ) for q ^ integers. 



(C13) 
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/co B t(l,0,g J r)--a 




FIG. 6. N umer ical evaluation of 1/r divergent part of / con t(l,0,g, r) f or r = 10 5 . Solid line repr esents the data obtained 
using Eq. (C13l, while the dots are obtained using the exact expression (C3l plugged in into Eq. (CI I. The plotted function is 



even in q and we only show the region q > 0. 



For positive integer q = N, one can not expand the function tp( — g(l — r)) from Eq. (CIO) at small r due to the 
divergence of i/^ - N). In Fig. \& we show a comparison between the exact result and its small r expanded form for 
the 1/r divergent part of 1(1, 0, q, r). 

We finally notice that Eq. ( |C10[ ) diverges, once p approaches an integer. Due to the parity of Eq. ( CIO ) with respect 
to p, we can restrict to the region p > 0. It is then convenient to rewrite i/j(p) + ip(—p) — 2tp(p) + ir cot(irp) + 1/p. As 
a result, the divergence arises due to the cotangent term. Therefore, the most divergent part of Eq. (CIO) at r — > 
and p — > N with N being a positive integer, is 



J(p,t) 

4tt 2 



8^2 - -^ncot{Trp) + R(p), 



(C14) 



where R(p) is the r-independent part in Eq. (CIO), contained in O(t ). Direct expansion of Eq. (C3| shows that it has 



a divergence of the type N /[4(p— N) 2 ] for any integer N, so R{p) can be rewritten as R{p) = Y^n=i m -ao 2 + • • • = 

^i\p->{Y — p) + . . ., taking into account only the most divergent terms 0[(q — 7V)~ 2 ], but not the terms that diverge 
as 0[(q — iV) -1 ]. Here . . . denotes other subleading divergent terms when q —J- N. Therefore, 



J(N + h,r) _(N + hf (N + hf 



At: 2 



8t 2 



4r 



0{h) 






[l + 0(h/N)} 



N 2 N 2 

..- [ l + { h/N) ] --[l 



■0{h/N)) + ^[l 



0(h/N)} 



(C15) 



Using Eqs. (CI) and (C15) we recognize h = —Nt, and therefore 

N N 2 + 2N + 2 
/ cont (l,0,iV,r) = 



8r 2 



+ G(r- 1 ). 



(C16) 



In Fig. [7|we show a numerical comparison between the obtained results. 

In order to get the C(r _1 ) term in the previous formula one should evaluate the 0(t 2 ) terms in (C7) and (C9) 
that we have not succeeded. However, we have verified the following two results 



leant (1,0,1,T) = 

J cont (l,0,2,r) 



5 

5 

: 4r2 



1 

4 



+ O(t ), 

+ O(t ). 



(C17) 
(C18) 
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8r 2 / cont (l,0,q,T) 

25 

20 

15 

10 

5 




FIG. 7. Numerical evaluation of 1/r 2 divergent part of I C ont(l, 0, q, r) for r = 10 -3 . For non-integer q values, 8r 2 / con t(l, 0, q, r) 
equals q 2 , while f or in teger q it is q 2 + 2q + 2, a s indicated by the dashed lines. The dots represent values obtained using the 
exact expression (C3l plugged in into Eq. (CI I at particular values q — 0.5,1,1.5,2,2.5,3. By solid line we show the result 



(C14l plugged in into Eq. (CI I. The inset shows zoomed result around q = 1, and we notice a very good agreement between 
the results | |C3[) and ( |C14[ ). For special q values such that p = q(l — r) approaches unity or any other integer, we notice the 
divergence of r 7 CO nt(l, 0, q, r). 



Appendix D: Integral I{r,a,q,r — 0) 



In this Appendix we consider the integral (41) from the main text at r = 0, i.e., 

'(x — r) 2 + a 2 y 2 



If 1 

I(r, a,q,r = 0) =— -r. d 2 xd 2 y ^ 

4tt 2 J [ a 2 + ( x -y) 2 ] 2 



(y-r) 



(Dl) 



and evaluate it for q — 1 and in the limit q — > 0. We notice the invariance /(r, a, q, 0) = I(r /a, 1, q, 0) and therefore 
we set a = 1, recovering a in the final result. 



1. The case I(r,a,q — 1, r = 0) 



We start with the integral (Dl| in the case q = 1. We represent the two-dimensional vectors x and y in polar 



coordinates. After performing the two angular integrations, and then the integration over x, the result can be expressed 
in the following form 



/(r, 1,1,0)= / dy[ J1 (r,y)+j 2 (y)+j 3 (r,y)], 
Jo 



where 



3i(r,y) 



1 



2y(y 2 + 4) 

h(y) = - 



2-y 2 

2 



(y 2 + l)(^ 2 -2r 2 -2) 
^/r 4 -2r 2 (y 2 -l) + (y 2 + l) 2 

-arctanh ' " 



J3(r,y) 



y 2 {y 2 +4) 3 / 2 



( y 2 + 4)3/2 \ y 2 + 2 

(y 2 + l)(r 2 y 2 + 2y 2 + 2r 2 + 2) ^ 
^ r 4 _ 2r 2 (2/ 2 _ !) + ( y 2 + ij2 



arctanh 



y V 7 ?/ 2 + 4 
y 2 + 2 



(D2) 



(D3) 

(D4) 
(D5) 
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After introducing a new variable t = y 

~ + l 



2, one easily performs the integration of ji(r, y) and obtains 



,00 j 

/ dyji(r,y) = - - 



In 



3r 2 + 9 



8v/(r 2 + l)(r 2 + 9) 



In 



(r 2 + 5)^(r 2 + l)(r 2 + 9)+ r 4 + 10r 2 + 17 



2 ln(r) + - In 2 - 



(9 



The function ,72(2/) is r- independent and gives a constant contribution 

dyhiv) = -In 2. 



(D6) 



(D7) 



The remaining yet unevaluated integral is L dyjj,{r : y). We notice that it contains a r- independent part [the term 
—2 from the square brackets in Eq. (D5|], which is necessary to make js(r,y) nonsingular at y — > 0. From the square 
root of the denominator of Eq. (D5 ), we conclude that ja(r, y) is sharply peaked when y takes values around r, with 
the height of the peak equal to ln(r) in the limit of large r. Therefore, we split j^ir, y) in the following form 



J3 (r, y) = J31 (r, y) + J32 (r,y), 



where 



r 2 arctanh 



3ai(r,y) 



y 2 +2 



y v / r 4 -2r 2 ( 2 / 2 -l) + (y 2 + l) 2 



(D8) 



(D9) 



In the last equation we keep y in the denominator in order make J32 (?*)?/) decaying faster than 1/y at large y. We 
emphasize here that the special choice (D9) for the first term of Eq. (D8) is not unique. For example, another 

2r 2 lnfa+1) 



possibility would be 



J / A /r*-2r2(y2-l) + ( !y 2 + l)2 



which is equally good, as it is integrable and describes well the function 



js(r,y) around its maximum. For simplicity, in the following we use the choice given by Eq. (D9). 

The subleading part of js(r, y) contained in J32(r,y) = js(r,y) — jai(r,y) can be safely expanded at larger, as it s 
main contribution in the integral comes from small values of y, i.e., from the region y <C r. Using Eqs. (D5|, (D8|, 
and ( D9 1 one obtains 



J32(r,y) = 



2/ 2 + 3 



1 



therefore 



(y 2 + 4) 3 / 2 y 



dy.J32(r,y) 



arctanh 



II 



vV + 4' 



o 



1 



7T~ 

~6~ 



In 2 



O 



(D10) 



(Dll) 



We perform the remaining integration L dy J3i(r,y) in the following way: We first calculate 



h\{r,y) 



dr 



J3i(r,y) 



arctanh 



iy/y 2 +i 
y 2 +2 



2y 



In 1 



(l + r 2 -y 2 + v / r 4_2r 2 (y 2 -l) + ( 2 / 2 + l) 2 ) 



(D12) 



Then, we expand J3i(r,y) at large r at lowest order in r and get 

arctanh 



Jli(r,y) = 



y\fy 2 +j 
y 2 +2 



In 1 



II 



(D13) 



where by . . . we denote the subleading terms at large r. The expansion of Eq. (D13) is a good approximation of 



J 31 (r,y) of Eq. (D12| only for y 



< 



For y > r, the function J3i(r,y) sharply drops to zero, contrary to its 



expanded form (D13), as one can see from the expansion of J3i(r, y). In the vicinity of r, at leading order one obtains 
J31 (r, y = r + 5) = (In r/r) [1 — 8/ In r] for 8 <C 1. Such an expansion determines a very large slope 0(1/ In r) for the 
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deviation of the function Jsi(r,y) around the point y = r. Therefore, we integrate the expanded result (D13l over y 
in the interval [0,r] and get 



dyJ3i{r,y) =-Li 2 



-r (r+ \Jr 2 + 4 J 



ln(r) — — ln(r) In 



+ ry/r 2 + 4 + 2 



7i~ 

~6 



= In 3 r + ^- In r + O 



(D14) 



The sublcading terms of Eq. (D14|, denoted by . . ., originate from the subleading terms of Eq. (D13) which do not 



change the stated result in Eq. (D14), as one can check, e.g., numerically. Finally, assuming the following form 

2 



Jo dy J31 (r, y) — A\n r + B In r + C + . . . , one obtains 



dr I dv — = — In r A In r + C In r ■ 

' y r 3 2 



(D15) 



After comparing the last expression with Eq. (D14) one gets A = 3,B = 0, and C — ir 2 /6. Combining this result 
with Eqs. (D8) and (Dll) one obtains 



f°° 2 

/ dyj 3 (r,y) =3 In r- 
Jo 



ln2 „(\ 



(D16) 



Collecting the obtained results (D6), (D7), and (D16l and using (D2|, one obtains 



I(r,a, 1,0) =31n 2 f-^) - 2 In / ' 



O 



(D17) 



where we have recovered the parameter a. 



2. The case I(r, a,q ->• 0, r = 0) 



The leading order term at q — > of Eq. (44) is given by I\{r) of Eq. (45). We could evaluate it using a procedure 



similar to the one employed during the evaluation of I(r, a, 1,0), see Appendix Dl First, one performs the two 
angular integrations, where one employs the following non-elementary integral 



/•2-7T 

/ dO ln(c - cos 6) = -2tt In ( 2c - 2 Vc 2 - 1 J , 



(D18) 



which can, e.g., be obtained by making use of the Jensen formulaP2l After performing one spatial integration, one 
obtains the result of the form 



h(r) = / dy 



y2 + 2 =arctanh[^±i)- 2 /ln( y 2 + l) 



V¥Tl 



y 2 + 2 



In 



V 2 + y 2 + 1 + v/r 4 - 2r 2 (y 2 - 1) + (y 2 + l) 2 ' 



2(y 2 + 1) 



(D19) 



For the last integral one could use a procedure similar to the one employed above for the integral (D2). The final 
result reads 



7 1 (r)=ln 2 (r)+ln(r) + O(r ) 



(D20) 



and therefore 



l(r,a,q,0) = q 2 



In 2 



©+'<) 



O 



0(q 4 ). 



(D21) 



Appendix E: Connection between the finite-a and 
the dimensional method 



the limit a —¥ is the same as f — > oo. The dimensional 



1. From poles to logarithms 



Let us define I(f,r) := f Ar I(r, a,q,r) and f = r/a 
[compare with the expression for W2 of Eq. (40)]. Hence 
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method gives, at fixed r > 0: 

'h(q) . h(q) 



J(r,T) 



~4t 

r 



0(1) 



(El) 



where 62 (g) and 61 (g) are even functions of g with a reg- 
ular Taylor expansion in q at g = 0. On the other hand, 
from the finite-a method, we know that 

I(r,T = 0) -^oo B 2 (q) In 2 f + B^q) lnf + 0(1)(E2) 

To match the two, we first observe that the integral 

(rd f ) 2 I(?,T) (E3) 

is not divergent and can be calculated at r = 0, giv- 
ing 2 -62 (g) from Eq. ( E2 1 , or directly calculated from 
Eq. ( |E1[ ) in the limit r — > is gives 1662(g)- Identifying 
the two coefficients, one gets 



B 2 (q) = 86 2 (g), 



(E4) 



which was used in the main text. 

When 62(g), -82(g) are non zero one cannot get in full 
generality a universal result for Bi(g), 61(g). It is easy to 
see since by simply changing the cutoff a by a finite scale 
-Bi(g) changes. However in the present case we can use 
the extra parameter, g, and we note that 62(g) = g 2 /8. 
Hence the integral 



I(r,r)={l--q 2 d 2 q )\ q=0 I(r 7 T) 



(E5) 

where 61(g) = (1 — i\q 2 d 2 ) | g= o6i(g) has a Taylor expan- 
sion starting at 0(g 4 ). In the finite-a method this new 
integral can have only a logarithmic divergence at large 
f, namely: 



I(f,r) =i?i(g)lnf + 0(1) 



(E6) 



and applying the same reasoning as above to the finite 
integral fdfl(f, r) we obtain 

S 1 (g)=46 1 (g). (E7) 

Hence in conclusion if 62(g) = g 2 /8 we conclude that: 



J(f,T = 0) 



g [In r + 7 (lnf)] 



(E8) 
(E9) 



+ 46i(g)lnf + 0(1), 
where 7 is non- universal, as used in the text. 

2. An illustrative example 



In this Appendix we study an example the connection 
between the two approaches used, the first one which 
keeps a finite ultraviolet cutoff a but uses r = 0, and the 
second dimensional method where a = but r > 0. In 



order to understand the connection, let us consider the 
following example, 



K(r, a,r) 



d z x- 



1 



\x\<l (x 2 + a 2 ) 



T37, (E10) 



which, although being simple, resembles the structure 
of the terms we encountered in Eqs. (40) and (41) and 



gives the essence of the diff erenc e of the two methods. 
The exact evaluation of Eq. (E10) is elementary and one 
obtains 



K(r, a, t) = 2-/T 



^ 2T (L 2 + a 2 y-a 2 
2t 



(Ell) 



Now we evaluate Eq. (E10) using the finite-a method. 



There, one first expands the integrand at small r, fol- 
lowed by the integration. Doing these two steps, one 
obtains 



K a (r,a,T) 



d 2 s 



1 



0(r) 



\x\<L 



L 2 



=27rln[ 1 + — ) +0(r). (E12) 



In the dimensional method, one first expands the inte- 
grand at small a, and then performs the integration. For 
r > 0, this procedure leads to a convergent integral, since 
the short-distance divergence is avoided by keeping r fi- 
nite. In such a way one obtains 



A- dim (r,a,r) 



1 



=2 



'©' 



\x\<L 

L 2 ^ 

~27 



d 2 a 
+ 0(a 2 ) 



— + 0(« 2 ) 
(E13) 



In the final step of the dimensional method one expands 
at small t the result (E13), finding 



K ailn (r,a,T) = 2n 



a 



1 
2t 



0(t) 



+ 0(a'- ZT ) 



(E14) 



The obvious discrepancy, at small a, between the two re- 



sults of Eqs. (E12) and (E14), is contained in the term 



tt/t in the latter; it arises due to setting a = in the 
integrand of Eq. (E13), which is reminiscent of setting a 
to zero in the last term of the exact result (Ell). As a 



consequence, the divergence at r — » in Eq. (E14) re- 



mains, which would have been canceled had we kept a 
finite a in the dimensional method. Therefore, in order 
to compare the two results of the finite-a and the dimen- 
sional method, one must neglect all the terms that are 
divergent in the limit r — > in the final result of the 
dimensional method. 

We note finally that expressions of the form a T , that 
appear in the integral (E10) and in similar ones, are 



treated differently in the finite-a method and in the di- 
mensional one. In the former one encounters the limit 
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t = and therefore a T = 1. In the latter one considers 
the limit a — > 0, and therefore a T = 0. This difference 
must taken into account if one want to compare the re- 
sults obtained by the two methods. 



Appendix F: From moments to distribution 



In principle the knowledge of the cumulants ( 67 ) of the 

9(r) 



ows to 



relative phase displacements 6 = 9(r) — 9(0) al i 
learn about the probability distribution function (PDF) 
of 0, which could be measured directly, e.g., in numerical 
simulations. 

The presence of poles in rj(q) at q — 1, and presumably 
integers, as well as the (related) mechanism of screening 
discussed in the introduction precluded us for now to 
have a complete knowledge of A(q) and 77(g) for all q, 
hence we are not able to fully characterize the PDF. 

In a more modest attempt, hopefully avoiding some of 
these problems, let us focus on the correlation at imagi- 
nary q = —ip, for which our result reads: 



(eP ) ~ exp ( p 2 t 2 £ 2 + [cp 2 - g(p)T 2 }£) 
g(p) = P 2 [^1e + i>(ip) + ip(-ip)], 



(Fl) 
(F2) 



where £ = ln(r/a) and g(p) is a positive even function, 
increasing for p > 0, which behaves as 0(p A ) at small p 
and as 2p 2 Inp at large p. Note th at th ere are no poles 
on the real p axis. At this stage (Fl) is just the gen- 



crating function of all cumulants of O and the explicit 
form (Fl) provides a resummation of the Taylor series 
in some vicinity of p = 0. It is not obvious whether 
this equality holds more globally, i.e., whether additional 
non-perturbative terms are also present, as is the case 
along real q. However it is maybe less likely for real p. 



To test that we must verify that (Fl ) is first of all an 
increasing function of p. Clearly, since g(p) ~ 2p 2 lnp 
this fails at fixed £ for some large enough p. So that for- 
mula cannot extend to r 2 lnp > cl + t 2 £ 2 . It is probable 
that for such large p one leaves the domain of validity of 
the (renormalized) perturbative expansion and a differ- 
ent calculation must be performed (such as an instanton 
calculation for real p) . 

Let us point out an interesting interpretation. Let us 
rewrite, to the same accuracy at large £: 

J^) ~ exp (cp 2 £ + p 2 T 2 £ 2 e- h ^l 1 ^ (F3) 

where h(p) — g(p)/p 2 ~ 0(p 2 ) at small p and h(p) ~ 
2 lnp at large p. From the large p> 1 behavior we can 
rewrite: 

(^«y ~ p>>1 exp (cp 2 £ + r 2 ^V~ f ) , (F4) 

i.e., indicating some deviations from the Gaussian behav- 
ior. 




FIG. 8. The correlation function C(q) for T = 1/5 (thick, 
green, dotted), T = 1/20 (red, thick, solid), and in the limit of 
T — > (black, thin, dashed). In the latter case it approaches 
C(<z)| T =o = 2sin(g/2)/g. 
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FIG. 9. The probability distribution JG6j), for T = 1/100 
(blue, dashed), T = 1/20 (red, solid), and T = 1/5 (green, 
dotted). 



Appendix G: A toy model 

Here we give results for a toy model which shows oscil- 
latory behavior in C(q) as a function of q. We write the 
energy for a particle in a parabola plus a cosine, which 
has a random shift £ € [—1/2, 1/2], 

H[9} = ^6 2 -g a cos(2Tr{0 + Z}) . (Gl) 

For simplicity we consider the limit of go —} 00. This 
restricts 9 to — £ plus an integer m. The expectation of 
e lg6 , given £ can then be written as 



<e^) 



E 



a iq(m—£)- 



(m-Of 



m— — 00 



E 



771 — — OO 



■^$ 3 (n((~ t qT),e- 2 * 2T ) 
7? 3 «,e- 2w2T ) 



(G2) 
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where ^3 is the elliptic 9 function and £ = — £. Then 



C(q) = (e^) = I d£ (e^ e ) 

1-1/2 



1/2 



(G3) 



We have plotted the result in Fig. [Hj For T — > 0, we can 
restrict the sum in Eq. (G2) to m = 0. This yields 



C(q) 



1 / 2 2 /ox 

d£e ige = -sm(") (G4) 

/-. ./_i/ 2 <Z V2^ V ; 



The corresponding probability distribution is 



P{6) 



e 



T=0 



\<*<\ 



(G5) 



This means that is uniformly distributed between —1/2 
and 1/2. For higher temperatures, the distribution will 
be smeared out. It then reads 



P(9) = 



e 2T 



This is plotted in Fig. [9] 



(G6) 



1 J. L. Cardy and S. Ostlund, Phys. Rev. B 25, 6899 (1982). 

2 J. Toner and D. P. DiVincenzo, Phys. Rev. B 41, 632 
(1990). 

3 T. Hwa and D. S. Fisher, Phys. Rev. Lett. 72, 2466 (1994). 

4 D. Carpentier and P. Le Doussal, Phys. Rev. B 55, 12128 
(1997). 

5 T. Giamarchi and P. Le Doussal, in Spin glasses and ran- 
dom fields, edited by A. Young (Singapore, 1997). 

6 M. P. A. Fisher, Phys. Rev. Lett. 62, 1415 (1989). 

7 G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. 
Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

8 T. Nattermann and S. Scheidl, Adv. Phys. 49, 607 (2000). 

9 D. Feldman, Int. J. Mod. Phys. B 15, 2945 (2001). 

10 M. Tissier and G. Tarjus, Phys. Rev. Lett. 96, 087202 
(2006). 

11 A. A. Fedorenko and F. Kiihnel, Phys. Rev. B 75, 174206 
(2007). 

12 L. Radzihovsky and J. Toner, Phys. Rev. B 60, 206 (1999). 

13 T. Giamarchi and P. Le Doussal, Phys. Rev. B 52, 1242 
(1995). 

14 P. Le Doussal, Int. J. Mod. Phys. 24, 3855 (2010). 

15 D. S. Fisher, Phys. Rev. Lett. 78, 1964 (1997). 

16 T. Klein, I. Joumard, S. Blanchard, J. Marcus, R. Cubitt, 
T. Giamarchi, and P. Le Doussal, Nature 413, 404 (2001). 

17 S. Bogner, T. Emig, and T. Nattermann, Phys. Rev. B 
63, 174501 (2001). 

18 A. Fedorenko, P. Le Doussal, and K. Wiese, (unpublished). 

19 Z. Ristivojevic, P. Le Doussal, and K. J. Wiese, Phys. 
Rev. B 86, 054201 (2012). 

20 A. Perret, Z. Ristivojevic, P. Le Doussal, G. Schehr, and 
K. J. Wiese, Phys. Rev. Lett. 109, 157205 (2012). 

21 Q. Zhang and L. Radzihovsky, Phys. Rev. E 87, 022509 
(2013). 

22 Q. Zhang and L. Radzihovsky, Europhys. Lett. 98, 56007 
(2012). 

23 P. M. Chaikin and T. M. Lubensky, Principles of condensed 
matter physics (Cambridge University Press, 1995). 



P. Nozieres, in Solids Far From Equilibrium, edited by 
C. Godreche (Cambridge University Press, Cambridge, 
1992). 

Note that since it involves order 2m in perturbation theory 
in the disorder, i.e., Aim ~ <? 2m , the result (3| holds only 
in the limit of r — > 00 at fixed q. For a large but fixed r, it 
will be more difficult to measure the behavior for larger q. 
M. Bauer and D. Bernard, Europhys. Lett. 33, 255 (1996). 
P. Le Doussal, K. Wiese, and Z. Ristivojevic, (unpub- 
lished) . 

S. Guruswamy, A. LeClair, and A. W. W. Ludwig, Nucl. 
Phys. B 583, 475 (2000). 

P. Le Doussal and G. Schehr, Phys. Rev. B 75, 184401 
(2007). 

A real crystal has a finite correlation length £b for the bulk 
translational order. Here we assume delta correlated dis- 
order and therefore capture physics at length scales larger 
than £b. 

T. Giamarchi, Quantum Physics in One Dimension 
(Clarendon press, Oxford, 2003). 

U. Schulz, J. Villain, E. Brezfn, and H. Orland, J. Stat. 
Phys 51, 1 (1988). 

P. Le Doussal, Ann. Phys. 325, 49 (2010). 
J. Zinn- Justin, Quantum Field Theory and Critical Phe- 
nomena (Clarendon Press, Oxford, 2002). 
P. Le Doussal, Europhys. Lett. 76, 457 (2006). 
V. Dotsenko and L. Fateev, Nucl. Phys. B 240, 312 (1984). 
V. Dotsenko and L. Fateev, Nucl. Phys. B 251, 691 (1985). 
V. S. Dotsenko, M. Picco, and P. Pujol, Nucl. Phys. B 
455, 701 (1995). 

R. Guida and N. Magnoli, Int. J. Mod. Phys. A 13, 1145 
(1998). 

A. Pmdnikov, Y. Brychkov, and O. Marichev, Integrals 
and Series: Volume 3: More Special Functions (Gordon 
and Breach, Amsterdam, 1998). 
http://functions.wolfram.com/07.27.17.0018.01. 
L. V. Ahlfors, Complex Analysis (McGraw-Hill, Inc., 
1979). 



